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t^ - , Abstract 

(N . 

Based on the R-parity violation option of the minimal supersymmetric Standard Model, we exa- 
mine the scenario where the massive gravitino, relic from the hot big-bang, is the lightest supersym- 
^ ' metric particle and can decay through one or several of the trilinear R-parity violating interactions. 

\& , We calculate the rates of the gravitino decay via the various three-body decay channels with final 

states involving three quarks and/or leptons. By taking into account the present constraints on the 
trilinear R-parity violating coupling constants and assuming the gravitino and scalar superpartner 
masses do not exceed O(10TeV), it turns out that the gravitinos could easily have decayed before 
the present epoch but not earlier than the big-bang nucleosynthesis one. Therefore, the considered 
scenario would upset the standard big-bang nucleosynthesis and we conclude that it does not seem 
to constitute a natural solution for the cosmological gravitino problem. 

PACS numbers: 04.65+e, 95.30.Cq, 98.80.Ft, 12.60. Jv, i4.8O.Ly 
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1 Introduction 



In supergravity theories ffl], the gravitino, namely the spin-3/2 supersymmetric partner 
of the graviton, weakly interacts with all the particle species (including itself) due to the 
small gravitational strength coupling \fG^ = 1/Mp, G^ being the gravitational constant 
and Mp the Planck scale. Hence, the gravitino-gravitino annihilation rate is extremely 
small so that the gravitinos should decouple at an early epoch of the universe history, and 
moreover at an epoch characterized typically by a temperature T higher than the gravitino 
mass: kT > m 3 / 2 [k being the Boltzmann constant) g. Therefore, the relic abundance of 
the gravitino should be large, which is often denoted in the literature as the cosmological 
"gravitino problem" . One of the first solutions for the gravitino problem to be envisaged is 
to compensate the large gravitino relic abundance by a gravitino mass sufficiently small, 
namely 777.3/2 < IkeV ||, to respect the limit on the present universe energy density: 
Q ~ 1. For heavier gravitinos, a second type of available solution is by shortening their 
lifetime so that they do not survive out at the late epochs 0. This can be realized in two 
characteristic options: either the gravitino is not the Lightest Supersymmetric Particle 
(LSP) and thus can decay into an odd number of superpartners through gravitational 
and gauge interactions (both of these ones couple an even number of superpartners), or 



it is the LSP. A gravitino LSP can be realized within various supersymmetric models, 
including the gauge mediated supersymmetry breaking models 0, H, the models of low 
fundamental energy scale |j| and even the conventional gravity mediated supersymmetry 
breaking models (for some specific set of the supersymmetry breaking parameters). If 
the gravitino is the LSP, it can decay only into the ordinary particles of the Standard 
Model. Such a decay channel must involve [[?]] both gravitational and the so-called Im- 
parity symmetry ||, || violating interactions (the latter ones couple an odd number of 
superpartners) . The R-parity violating {ftp) interactions are written in the following 
superpotential, in terms of the left-handed superfields for the leptons (L), quarks (Q) and 
Higgs of hypercharge 1/2 (H) and the right-handed superfields for the charged leptons 
(E c ), up and down type quarks (U C ,D C ), 

W fl P = £ ^jkUL 3 El + X' ljk UQ 3 Dl + l -\'l ]k UtDpi + fjuHLi) , (1) 

i,j,k being flavor indices, ^ijk, Kjki Kjk dimensionless coupling constants and fii dimen- 
sion one parameters. Note that in the scenario in which the gravitino is not the LSP, the 
gravitino preferentially decays into an ordinary Standard Model particle and its super- 
partner through gravitational interactions, since the ft p coupling constants are severely 



constrained by the low-energy experimental bounds obtained at colliders (10], [II]. 

In the scenario of an unstable gravitino heavier than the LSP, the gravitino decay 
can produce an unacceptable amount of LSP which conflicts with the observations of 
the present mass density of the universe JE3] . The scenario containing an unstable LSP 
gravitino, having decay channels which involve ftp coupling constants, is based on the 
violation of the R-parity symmetry. Now, neither the grand unified theories, the string 
theories nor the study of the discrete gauge symmetries give a strong theoretical argument 
in favor of the conservation of the R-parity symmetry in the supersymmetric extension of 



the Standard Model fllP] . Hence, the scenario with an unstable LSP gravitino constitutes 
an attractive possibility which must be considered as an original potential solution with 
respect to the cosmological gravitino problem. 

Our main purpose in the present work is to determine whether the scenario of an 
unstable LSP gravitino decaying via ftp interactions constitutes effectively a natural 
solution to the cosmological gravitino problem. We will concentrate on the trilinear 
ftp interactions, namely X^LiLjE^ \' i ^ k LiQjD'f, and X'l^UfD^D^. With the goal of en- 
hancing the gravitino instability, we will consider an optimistic type of scenario in which 
several trilinear fl p coupling constants have simultaneously non-vanishing values. We will 
not consider the bilinear ft v term faHLi of Eq.([l|) which can always be rotated away 
by a suitable redefinition of superfields ]nj. The bilinear ftp interactions as well as the 
possible alternative of a spontaneous breaking of the R-parity symmetry have been con- 



sidered, within the context of the cosmological gravitino problem, in a recent study [11] 
which examines the two-body gravitino decay mode into photon and neutrino (we will 
compare the results of [14]] with ours). We also mention here the works [[TJ|] and |16|] in 



which have been considered the ftp decay reactions of single nucleon into a gravitino and 
a strange meson. 

In Section §, we determine the constraints arising from the cosmological gravitino pro- 
blem in a scenario characterized by an unstable gravitino. In Section |||, we calculate 
the dominant gravitino decay rates involving the trilinear ft v interactions. Finally, in 
Section |, we give a quantitative discussion aimed at determining whether the scenario 
of an unstable LSP gravitino decaying through ftp interactions can effectively constitute 



a natural solution with regard to the cosmological gravitino problem. This quantitative 
discussion is based on the results of Sections and |3| as well as the present low-energy 



constraints on the fL, coupling constants [II 



2 Constraints on the gravitino decay temperature 

In this section, we derive the constraints |J characteristic of a scenario containing an 
unstable gravitino supposed to solve the cosmological gravitino problem. 

At the temperature kT « 7713/2, the expansion rate of the universe is H(kT « 7/13/2) ~ 
y/ m 3/2(kT) 3 / 2 jMp = m\, 2 /Mp. Besides, the gravitino decay having the larger rate is the 
two-body decay into a particle of the Standard Model and its supersymmetric partner. 
The largest gravitino decay rate is thus of order, 

Tmax ~ ml /2 /M P , (2) 

which is smaller than H(kT « 777.3/2) by a factor Mp/ 1713/2- Therefore, the equilibrium 
condition, 

r > H(T), (3) 

T being the gravitino decay rate, can only be reached for temperatures much smaller than 
kT m m 3 / 2 . In other words, the gravitinos cannot decay as long as the temperature has not 
dropped far below the gravitino mass. This conclusion has two important consequences. 

The first consequence is that the gravitino decays occur at temperatures at which the 
gravitinos cannot be created in collisions, so that the gravitino population can effectively 
be reduced by the gravitino decays. 

The other consequence is that the equilibrium condition (Eq.(|3|)) is reached at an 
epoch at which the cosmic energy density is dominated by the gravitino energy density 
P3/2. Therefore, in determining the temperature T3/2 at which the equilibrium condition of 
Eq.(^|) is reached (the gravitino decay temperature) by solving the equation T = H(T 3 / 2 ), 
one must take the following expression for the expansion rate, 



H(T) M f-^p m . (4) 

Since the gravitino mass density reads as, 

3C(3) g(T) 3 

P3/2 = —^r^f^ m ^ kT ) ' ( 5 ) 

where ((x) is the Riemann zeta function and £(3) = 1.20206 . . ., the equation Eq.@ can 
be rewritten as, 



H(T) 



mm^E m sn (6) 

\ 7T g[T d ) M P 



In Eq.(||) and Eq.(f|), g(T) counts the total number of effectively massless degrees of 
freedom (those of species with mass m <C kT), and Td is the gravitino decoupling tem- 
perature, namely the temperature at which is reached the out of equilibrium condition 
T an < H(T), T an being the gravitino-gravitino annihilation rate. The gravitino decou- 
pling temperature is typically larger than the gravitino mass, as we have mentioned in 



the introduction. The presence of the factor g(T)/g(Td) in Eq.fl5|) and Eq.(||) allows to 
take into account the increase of the temperature due to the conservation of the entropy 
per comoving volume S = (2ir 2 /4:5)g(T)(kT) 3 R 3 , where R = R(t) is the time dependent 
scale factor. From Eq.(^) it follows that the gravitino decay temperature is given by, 

kT ^ w (m)Ms) (^f) ■ (7) 

At the temperature T 3 / 2 , most of the gravitinos decay. After the gravitinos have decayed 
and their decay energy has been thermalized, the temperature rises to the value kTL 2 
given by, 

3/2 * V g(T 3/2 ) J K 

The equation Eq. @ is nothing else but the expression of the energy density conservation. 
From Eq.(0) and Eq.@ it follows that kTy 2 reads as, 

/45 V/- 1 g(T d ) 



kT ^{^) M^^^- <9) 

In order to avoid large relic abundance of an unstable gravitino, this one must decay 
before the present epoch. This requirement imposes the following bound on the increased 
gravitino decay temperature TL 2 , 

% j2 > 2.75K (kT^ /2 > 2.36 l(T 10 MeV). (10) 

Furthermore, the decay of the gravitinos leads to an increase of the temperature by a 
factor which can be deduced from Eq.flT]) and Eq.@: 

n /2 „ 45VV2C(3)V7 (T n^\- 1/9 f m m v /3 

1^ [g(T 3/ 2)^g(T d )) [—-) . (li) 
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In order to derive the bound of Eq. (|TTD , we have used the fact that T < T max (see 
Eq.(||)). This increase of the temperature leads to an increase of the entropy density 
s = (2ir 2 /4:5)g(T)(kT) 3 and hence to a decrease of the baryon-to-entropy ratio (or baryon 
number in a comoving volume) B = ub/s, ub being the baryon number density. There- 
fore, if the gravitinos decay after the nucleosynthesis epoch, the baryon-to-entropy ratio 
during the nucleosynthesis epoch must be much greater than the present one, leading 
to the production through nucleosynthesis of too much helium and too little deuterium 
(compared with the constraints on the primordial abundances derived from observational 
data) 0. In conclusion, if the gravitinos decay before the present epoch, those decays 
have to occur before the nucleosynthesis epoch, or in other words, the increased gravitino 
decay temperature Ty 2 has to be higher than the nucleosynthesis temperature, namely, 

kT^ /2 > OAMeV. (12) 

In summary, in a scenario of unstable gravitino protected against large gravitino 
relic abundance, the gravitino must decay before the present epoch, namely kTL 2 > 



2.36 10~ 10 MeV (see Eg. ([TDD), or even at an earlier epoch than the nucleosynthesis one, 
namely kTL 2 > OAMeV (see Eg. ([12])), if one reguires that this scenario does not conflict 
with the current cosmological ideas on nucleosynthesis. Hence, in order to determine 
whether a given scenario of unstable gravitino can effectively constitute a solution to the 
gravitino problem, one has to calculate, within this scenario, the gravitino decay rate. 
Indeed, the value of the gravitino decay rate gives, through Eg.(^), the gravitino decay 
temperature TL 2 , and thus allows to verify whether the constraints of Eq. fllOD and Eg. ([12]) 
are respected or not. The calculation of the gravitino decay rate will be done in the next 
section, within the scenario we consider. 

3 Rates of the gravitino decay via trilinear R-parity violating 
interactions 

As we have discussed in the introduction, in the scenario we study, the gravitino is the LSP, 
so that it can decay only via Rp interactions, and the R-parity symmetry is violated only 
by the trilinear terms of the superpotential |l|. Hence, in our scenario, the gravitino can 
only decay through the trilinear R p couplings of Eg. ([I]). The gravitino decay processes 
which involve the trilinear R p interactions of Eg.(|l|) and have the dominant rates are 
obviously the decays through a virtual scalar superpartner into three ordinary fermions 
involving one gravitational and one R p coupling. Those processes are represented in Fig.[I], 
Fig.g and Fig.|| for the R p interactions of type Xi jk LiLjE k , X' ijk LiQjDl and \'{ jk U^DjDl 
(see Eg. (|T])), respectively. The Feynman diagrams of Fig .[I], FigfJ and Fig.|3] have been 
derived from the structures of the relevant interaction lagrangians. Let us collect these 
relevant interaction lagrangians. First, the locally supersymmetric lagrangian term which 
describes the coupling between the fields of the gravitino {^p), an ordinary fermion (ip) 
and its scalar superpartner (0) reads as, 

C = "Vi" Wr9 ^ V + h-c '' (13) 

the L/R indices standing for Left/Right chirality, 7 M being the Dirac matrices and M* = 
(StiGn)" 1 ^ 2 = 2.4 10 18 Ge^ the reduced Planck mass. With our conventions, the reduced 
Planck mass M* and the gravitino mass 7723/2 are related through the formula, 

m ^ = tw: (14) 

where F denotes the scale of spontaneous supersymmetry breaking. Secondly, the relevant 
lagrangian of the R p interactions of type A, A' and A", originating from the superpotential 
D, is the following, 



'./'- 



C R P =J2~ X ijk - (piLekRe jL + e jL e kR u iL + e* kR ^ R e jL - (z *-> j) J 

~\jk [ ViLdkR,djL + djidkRViL + dlftViftdjL 

— ZiLdkRUjL — UjLdkRen — dl R e^ R UjL j 

+ h.c, (15) 



1 

ijk 9 I Lb iR u 'jR u 'kL ~l~ ^jRU'iR^kL 



-Kik 7T u*d ]R d c kL + 2d* R u iR d c 
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Figure 1: Feynman diagrams of the gravitino three-body decay processes involving the XijkLiLjE^, in- 
teractions. These interactions are represented by the circled vertex. The other vertex correspond to 
gravitational couplings. The plain (dashed) lines symbolize the propagation of fermionic (scalar) par- 
ticles. The arrows denote the flow of momentum of the associated particles. Finally, e, e, v, v and G 
represent respectively the charged leptons, the charged sleptons, the neutrinos, the sneutrinos and the 
gravitino, i,j, k are generation indices, the L 7 R indices stand for Left, Right (chirality) and the c expo- 
nent indicates a charge conjugated particle. We have not drawn the charge conjugated reactions which 
can be trivially deduced from these ones. 



where e, z/, d and u represent respectively the fields of the charged leptons, the neutrinos, 
the down and the up quarks, a ~ indicates the scalar superpartner and our notation is 



such that one has, for instance, v? R 



WM- 



As shown by Fig.[l], when the gravitino decays via a A^ coupling, it decays into two 



x. a 



charged leptons and a neutrino as G -^ UiejEk^iSjek- Based on the lagrangians of Eq. flT3[) 
and Eq.([HJ), we have calculated the rate of this gravitino decay reaction, which receives 
three contributions as can be observed in Fig.[L]. The full analytical result for the spin 
summed amplitude of this gravitino decay reaction is given in Appendix [FJ. We give here 

the analytical result for the integrated rate of the gravitino decay G -^ ViCjCk-, under the 
assumption that the masses of the final state particles are negligible with respect to the 
gravitino and exchanged scalar superpartner masses (we will comment on this hypothesis 
below), 
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where the mass functions Q(m), S(r7i 1 ,m 2 ), A(mi, 777.2) and Y,(mi,m 2 ) are defined in 
Appendix [A|, Sp(x) = Polylog{2,x) = Li 2 (x) is the dilogarithm or Spence function and 
m„ iL , mg. L and m^ kR are respectively the masses of the superpartners v iL , e^ L and e^R- 
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Figure 2: Feynman diagrams of the gravitino three-body decay processes involving the X' i j k LiQjD k 
interactions. These interactions are represented by the circled vertex. The other vertex correspond to 
gravitational couplings. The plain (dashed) lines symbolize the propagation of fcrmionic (scalar) particles. 
The arrows denote the flow of momentum of the associated particles. Finally e, e, v 1 i>, d, d, u, u and 
G represent respectively the charged leptons, the charged sleptons, the neutrinos, the sneutrinos, the 
down quarks, the down squarks, the up quarks, the up squarks and the gravitino, i,j, k are generation 
indices, the L, R indices stand for Left, Right (chirality) and the c exponent indicates a charge conjugated 
particle. We have not drawn the charge conjugated reactions which can be trivially deduced from these 
ones. 
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Figure 3: Feynman diagrams of the gravitino three-body decay processes involving the \'!- k UfD^D'{. 
interactions. These interactions are represented by the circled vertex. The other vertex correspond to 
gravitational couplings. The plain (dashed) lines symbolize the propagation of fermionic (scalar) particles. 
The arrows denote the flow of momentum of the associated particles. Finally, d, d, u, u and G represent 
respectively the down quarks, the down squarks, the up quarks, the up squarks and the gravitino, i,j, k 
are generation indices, the L, R indices stand for Left, Right (chirality) and the c exponent indicates a 
charge conjugated particle. We have not drawn the charge conjugated reactions which can be trivially 
deduced from these ones. 



The hypothesis that the masses of the final state particles, in the reaction G -^ ^je^, 
are negligible with respect to the gravitino and exchanged scalar superpartner masses 
has been made in order to obtain a more simple form for the analytical result of the 
integrated gravitino decay rate (see Eq. flTBD ). Nevertheless, in Section f|, the numerical 
results will take into account the full mass effect of the final state particles. Moreover, 
this assumption that the masses of the final state particles are negligible with respect to 
the gravitino and exchanged scalar superpartner masses represents a good approximation 
justified by the two following points. First, the present experimental limits on the scalar 
superpartner masses (noted generically by m), obtained at colliders within an ft p model 
containing a non-vanishing Yukawa coupling of type A, A' or A", are typically of order 



m ~ lOOGeF HITf . Secondly, as we will see, in the quantitative discussion of Section 



we will not need to consider gravitino masses smaller than 50GeV . 

Finally, we note that the integrated rate of the gravitino decay G -^ ViCjCk is equal to 



x. a 



i j k 



the integrated rate of the charge conjugated decay G — > z^e,efc, which is given by the 
formula of Eq. ([16]) . 



_ A 



The two kinds of gravitino decay channel involving A^- fe couplings are G -^ Uidjdk, Vidjdk 

and G — > eiUjdk, eiUjd^ (see Fig.|2|). Since the structures of the X^k and A^- fe interaction 
lagrangians are similar (see Eq.(|TJ), Fig .[I] and Fig.0), the integrated rate of the gravitino 

~ Kk - ~ Kk - r-?i 

decay G -4 v^djd^, (G -^ CiUjdk) is obtained from the formula of Eq.(nJ|) simply by 
changing \j k to A^- fc , adding an extra color factor of N c = 3 and replacing the masses 
as follows: m^ L -> m 9iL (rn SiL ), m SjL -> m SjL (m % J and m SkR -»■ m ikR {m SkB ). The 
approximation of neglecting the masses of the final state particles is also well motivated 
in the calculation of the gravitino decay rate involving \'^ k coupling constants, except 

~ A' — _ 

when the gravitino decays into the top quark as G -^ e^,d^e^d k . However, recall that 
in Section f|, the mass effect of the final state particles will be included in the numerical 
computations. 

In case the gravitino decays through a coupling of type A" jfc , it decays into three quarks 

A" 

as G -^ Uidjdk,Uidjdk, as illustrates Fig.[3]. In Appendix P, we give the full analytical 
result for the spin summed amplitude of this gravitino decay process. In the following, we 

present the analytical result for the integrated rate of the gravitino decay G -^ Uidjd kl 
in the hypothesis the masses of the final state particles are negligible compared with the 
gravitino and exchanged scalar superpartner masses (see above for detailed comments on 
this approximation), 
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[s P {—^) - s P (-^) + s P ( 3/2 2 "'« ) + Sp( 3/2 2 "-* ■ 

m ln. m \ R m l R m \ R 

(17) 

where N c = 3 is the number of colors and Q(m), E(mi, 777.2), A(mi, 777.2) an d S(mi,m 2 ) 
are defined in Appendix [A[ 



4 Numerical results and discussion 

In this section, we present the numerical results for the gravitino decay temperature TL 2 , 
within the scenario characterized by an LSP gravitino having a decay channel which 
involves trilinear ftp interactions. Then we will determine whether the values found for 
the temperature TL 2 respect or not the limits given in Eq. fllOD and Eq. ([l2|) . We will thus 
be able to conclude on the possibility for the scenario with an unstable LSP gravitino, 
decaying via A, A' or A" couplings, to solve the cosmo logical gravitino problem. 

• \jk couplings 

We first consider the case where one fl p coupling constant of type A^ is dominant so 
that all the other fi p coupling constants are negligible in comparison. The assumption 
of a single dominant fl p coupling constant is often adopted in the literature (as we will 
remark below) for simplification reasons. This assumption find its justification in the 
analogy between the structures of the fl p and Yukawa couplings, the later ones exhibiting 
a strong hierarchy in flavor space. In this first case, the gravitino decay temperature TL 2 
is derived (through Eq.(^)) from the gravitino decay rate involving A^-fc coupling constants 

Y{G -+ UiejCk), which is given in Eq.(|TB|). 

In Fig||(a) and Fig.[|(b), we present the gravitino decay temperature T^, 2 as a func- 
tion of the scalar superpartner mass fh for two gravitino masses, in the case of a single 
dominant fl p coupling constant of type A^. Let us detail some particular points con- 
cerning the derivation of those numerical results. First, for simplification reason, we 
have set all the scalar superpartner masses to a common value, denoted by fh. In or- 
der to maximize the gravitino decay rate, and thus the temperature Ty 2 (see Eq.(|])), 
we have also assumed that the single dominant fl p coupling constant is A233, which is 
one of the A^ coupling constants having the weakest present low-energy constraint, 
namely A233 < 0.06(m fR /100GeV) [0, [18]], and that A233 is equal to its present low- 



energy bound. We note that the present low-energy constraint on A233 has been ob- 
tained under the hypothesis that A233 is the single dominant ftp coupling constant 



which is consistent with our assumptions. At this stage, a particularity due the struc- 
ture of the Xijk couplings must be described. Since the A^ coupling constants are anti- 
symmetric in the i and j family indices, in the considered case where A233 is a dom- 
inant ftp coupling constant, A323 is a second dominant ftp coupling constant (having 

the opposite value) so that the main gravitino decay channels are G ^ u^rf, v^jt 

and G ^ v T nf, v T Jir. This point has been taken into account in the calculation. 
The decay temperature in Eq.(^) is proportional to the combination of statistical fac- 
tors g(Td) 1 / 4 / 'g(T 3 / 2 ) 1//2 • Let us here assume tentatively that the thermalized degrees of 
fredom at the gravitino decoupling and decay temperatures are the same as those for 
the minimal supersymmetric Standard Model and the present epochs, respectively. Li- 
sing then the values g(T d ) ~ 915/4 ~ 228.75, g{T 3/2 ) = 43/11 = 3.909, one obtains: 
Tg/ 2 oc g(Td) 1 / 4 / ' g{T^/2) 1 ^ 2 — 1.96. The actual numerical value of this factor might be 
larger but not by very much. In the numerical results, the above ratio of statistical fac- 
tors has been set to unity (we will come back to this point later). Finally, the gravitino 
decay rate has been multiplied by a factor of 2 in order to count the charge conjugated 
gravitino decay process and the masses of the final state particles have been taken into 
account in the computation. 
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Figure 4: Gravitino decay temperature kT^, 2 (in MeV) as a function of the scalar superpartner mass 
m (in GeV/c 2 ), with m 3 / 2 = 50GeV/c 2 in Figures (a,c,e) and m 3 / 2 = lTeV^/c 2 in Figures (b,d,f). The 
temperature T' /2 represented by a plain line in Figures (a,b), (c,d) and (e,f) has been obtained (via 
Eq.(||)) from the gravitino decay rates involving the ft p coupling constants A233 = 0.06(m/100Gey), 
A 331 = 0.48 for m = lOOGeV (see |H| for the numerical dependence on in) and A213 = 1-25, respectively. 
The temperature Tg, 2 represented as a dashed line in Figure (f) has been derived from the sum of 
gravitino decay rates involving the following ftp coupling constants, A' 132 = 0.34 for m = lOOGeV^ (see 
(ID for the numerical dependence on m), A 211 = 0.09(fh/l00GeV), A 223 = 0.18(m/100GeF), A 311 = 
0.1Q(m/100GeV), A121 = 0.05(m/100Gey) and A 233 = 0.06(m/100GeF). 

We observe in Fig.|](a) and Fig|4|(b) a decrease of the gravitino decay temperature 
T3 , 2 when the scalar superpartner becomes heavier (particularly in Fig.f|(a) for m 3 / 2 ~ 
50GeV). This is due to the fact that TL 2 is proportional (see Eq.(|])) to the square root 

of the sum of gravitino decay rates Y{G ^ v^tt) + T(G -4- v T [if) which decreases as the 
scalar superpartner mass increases. Nevertheless, this decrease of TL 2 when fh increases 

is attenuated by our choice of taking A233 = 0.06(m/100GeV), since T(G -^? v^tt) oc A233 
and T(G -^ v T {if) oc A| 23 = A233 (see Eq.(|TB])) so that TL 2 oc A233. We also remark by 
comparing Fig.^(a) and Fig.|](b) that TL 2 is getting higher if 7713/2 is enhanced. Once 
more, the explanation is related to the behavior of the sum of gravitino decay rates 

T(G -^> 3 v^tt) + T(G -^ z/ T /if): those ones of course increase with m 3 / 2 . We conclude 
from Fig.[|(a) and Fig.|](b) that the gravitino decay temperature TL 2 fulfills the condition 
Hi] but not 0, in case one fl p coupling constant of type \jk is dominant. This conclusion 



cannot be modified by the value of the factor g(Td) 1 / 4 / 'giT^i^) 1 ^ 2 of Eq.(|^) which is not 
likely to be significantly large. We must also note that only scalar superpartner masses 
smaller than 0(TeV) have been considered in Fig.|](a) and Fig.|](b). This is to ensure that 
the "hierarchy problem" , namely the problem of natural coexistence of the electroweak 
symmetry breaking scale and the scale of new physics underlying the Standard Model 
(grand unification scale, string scale . . . ), can be solved by the presence of supersymmetry. 
Now, if the gravitino mass (which is also the LSP mass), and thus the scalar superpartner 
masses, reach values as large as ~ 10 2 TeV, the condition of Eq.flT2|) can be fulfilled, in 
case one ftp coupling constant of type \j k is dominant. 

• Kjk couplings 

Secondly, we suppose that one ftp coupling constant of type X'^ k is dominant compared 
with all the other ft p coupling constants. In this framework, the temperature TL 2 is then 
obtained (through Eq.@) from the whole gravitino decay rate, which is the sum of the 

two gravitino decay rates involving a A^ fc coupling constant, namely T(G -^ Uidjdk) and 

T(G -4 eiUjdk) (see Section y). 

In Fig.[|(c) and Fig.f|(d), we show the temperature TL 2 as a function of the scalar super- 
partner mass m for two gravitino masses, in the case of a single dominant fl p coupling 
constant of type A^ fc . The explanations on the calculation of TL 2 concerning the numeri- 
cal factors and the masses of the particles involved in the gravitino decay process, which 
are given above for the case of a dominant A^ coupling, still hold in the present case. In 
order to optimize the temperature TL 2l we have also supposed here that the single domi- 
nant ftp coupling constant is A 331 , which is one of the \'^ k coupling constants having the 
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weakest present low-energy constraint, namely A 331 < 0.48 for rriq = lOOGeV |TI| (see |19j 
for the numerical dependence of this limit on the squark mass) , and that A' 331 is equal to 
its present maximum allowed value. The present low-energy constraint on A' 331 has been 

derived by assuming that A' 331 is the single dominant ftp coupling constant |BJ, which 

- \'. — 

is once more consistent with our assumptions. The gravitino decay rate T(G -^ rid) 



depends on the top quark mass that we have fixed at m t0 p = 174.2GeV p0| . 
The dependence of the temperature TL 2 on in and m 3 / 2 observed in Fig||(c) and Fig.f|(d) 
is similar to that seen in Fig.f|(a) and Fig.|](b), which has been described above. 
It appears clearly in Fig.|](c) and Fig.f|(d) that the gravitino decay temperature T 3 , 2 res- 
pects the constraint of Eq. fllOD but not of Eq. (|l2|) , in case one ftp coupling constant of 
type A^- fc is dominant. It must be mentioned at this level that, as before, scalar superpart- 
ner masses larger than 0(TeV) have not been considered in order to preserve the natural 
coexistence of electroweak symmetry breaking and new physics scales, guaranteed by su- 
persymmetric theories. We have found that for gravitino masses greater than ~ 10 2 TeV, 



the constraint [12] can be respected in the case of a dominant \'^ k coupling. 
• Kjk couplings. 
For the case of a single dominant ftp coupling constant of type A" fc , the temperature 



~ A" 



TL 2 is deduced from the gravitino decay rate involving X'^ k coupling constants T(G 
Uidjd k ) (see Eq.flT7|)). 

In Fig. |(e) and Fig.|(f) (plain line), we have represented TL 2 as a function of ih for two 
values of m 3 / 2 , in the case of a single dominant A" fc coupling constant. The explanations 
given above on the computation of T 3 , 2 concerning the numerical factors, the particle 
masses and the anti-symmetry of the A^ coupling constants, still applies in the present 
case since the A" k coupling constants are anti-symmetric in the j and k generation indices. 
We have also supposed here that the single dominant ftp coupling constant is A 213 , which 
is one of the X"j k coupling constants having the weakest present bound, namely A 213 < 1.25 
for a supersymmetry breaking scale of Msusy ~ ITeV |21|, |22|], and that A 213 is equal 
to its present limit. Let us make a few remarks on the present bound on A 213 . This 
constraint is not, as the previous ones, a low-energy experimental bound obtained at 
colliders: it originates from the requirement that A 213 remains in the perturbative domain 
(A 2 ^ 3 /47r < 1) up to the scale of gauge group unification. We also mention that this 
constraint depends weakly on the value of Msusy, and that it has been obtained under 
the hypothesis that A 213 is the dominant ftp coupling constant. 
We see in Fig.|](e) and Fig^(f ) (plain line) that the temperature T 3 , 2 fulfills the condition 



H] but not [LJ 

In Fig.[5] (plain line), we show the values of the gravitino and superpartner masses for 
which hTL 2 > OAMeV. These values have been obtained by assuming that the dominant 
ftp coupling constant is A 213 and is equal to its present limit. We see in Fig.|^ (plain 
line) that as the superpartner mass increases, larger gravitino masses are needed to have 
kTL 2 > OAMeV. The reasons are that T 3 / 2 increases with the gravitino mass but is 
suppressed if the superpartner is getting heavier (see previous figures). We can however 
conclude from Fig.[5] (plain line) that, typically, gravitino (and scalar superpartner) masses 
larger than ~ 80TeV are needed so that the condition of Eq. (|i~2|) can be fulfilled, in the 
case where A 213 is the dominant coupling constant. 

Finally, we consider the situation in which a maximum number of ftp coupling con- 
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Figure 5: Domains of the m 3 / 2 (GeV/ c 2 )-m(GeV/ c 2 ) plane (gravitino versus superpartner mass) in which 
the gravitino decay temperature is higher than the nucleosynthesis one, namely kT^, 2 > OAMeV. 
The region situated above the plain (dashed) line corresponds to kT^< 2 > OAMeV in case the do- 
minant fi p coupling constant is A 213 = 1.25 (the dominant fi p coupling constants are A' 132 = 0.34 
for m = lOOGeV, A' 211 = 0.09(m/100GeV), A 223 = 0.18(m/100GeV), A 311 = 0.10(m/100GeF), 
A121 = 0.05(rh/100GeV) and A233 = 0.06(m/100GeV r )). We note that the considered Jt p coupling cons- 
tants are fixed to the value of the perturbative limit (V47r) in the mass domains where their present 
bound exceeds this limit. Finally, the colored region corresponds to the situation m 3 / 2 > m which must 
be considered within a scenario where the gravitino is not the LSP. 



stants having the weakest present bounds are simultaneously dominant (and the consid- 
ered fl p coupling constants are equal to their present limit, as before). Based on the 
strongest constraints on the products of fl p coupling constants and on the review [11] 
of the present limits on the single fl p coupling constants, we have found that this opti- 
mistic situation corresponds to the case where the simultaneously dominant fl p coupling 
constants are, 

■^132' ^211; -^223' ^311' Al21 an d A 233 , 



which have the present bounds A' 132 < 0.34 for rrtq 



lOOGeV^ (and for instance A' 132 ~ 1.2 



for m- q = ITeV) 0, 0, A' 211 < 0.09(m^ /lOOGeV) pf, A 223 < 0.l8(m- b J100GeV) 



23j, A 311 < 0.10(m d - /lOOGeV) [0, g3|, Am < 0.05(m e - fl /100Ge1/) |10|, |f] and A 233 < 
0.06(m^ i j/100Gel / ) |T^, [1^]. At this stage, an important point to note is that the X^ k and 
\"j k couplings cannot be simultaneously dominant, since the experimental constraints on 
the proton decay rate force any product AL fc A^v fc / to be smaller than 1CT 9 , in a conser- 
vative way and for squark masses below ITeV. This result has been obtained in [|26| by 
calculating the proton decay rate at one loop level. In contrast, we have checked here that 
no strong constraint exists on any product taken among the ftp coupling constants A' 132 , 
A 211 , A 223 , A 311 , A121 and A 2 33. We also note that in the optimistic situation described 
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above, no other ftp coupling constant of type \' 22 k (-Msifc) can De supposed simultaneously 
dominant, since the present constraint on A' 223 (A 311 ) has been obtained under the hy- 
pothesis that A 223 (A311) is the only dominant coupling constant among the \' 22k (A 31fc ) 
coupling constants P3j] . This interpretation corresponds admittedly to a strong version of 
the single coupling constant dominance hypothesis implying that the contributions from 
the subdominant coupling constants do not interfere destructively in the observable used 
to infer the bound. Similarly, the present bound on A' 132 |~l], HH ( a coupling constant 
of type X' 32k [|K], |24|) has been derived by assuming that it dominates over all the A' ljfc 
(X' ilk and \' i2k ) coupling constants, so that no additional fi p coupling constant of type 
^ijk (^32/fc) can be assumed simultaneously dominant. Once again, if we use the present 
bound on A 211 , all the other X' 21k and X' uk coupling constants must neglected [jTSj , so that 



no additional ftp coupling constant of type X' 2lk can be considered simultaneously. A' 231 
also cannot be simultaneously dominant, for another reason which is the existence of the 
following constraint, 

A^A'^ < 4.7 10- 5 {m-JimGeVf g5f. 

Finally, the remaining A' coupling constants, X 23k {k = 2,3) and X 33k {k = 1,2,3), can- 
not be additional dominating couplings, because the present bounds on those have been 
obtained by assuming the dominance over all the A 2 „ fe and X' 3 j k couplings, respectively 
111 , HI]- Concerning the A^ coupling constants, the constraint on Ai 2 i (A 233 ) has been 



obtained under the assumption that A i2 i (A 233 ) dominates among the A^ and A' llfc (A^ 



and A 23 fc) coupling constants | 18fl . Therefore, in the optimistic framework described above, 
while the hypothesis that A121 and A233 simultaneously dominate remains consistent with 
the assumptions adopted in fl8|], no additional fl p coupling constant of type Xijk can be 



considered simultaneously. 

In Fig.|](f) (dashed line), we have represented T 3 , 2 as a function of m, in the optimistic 

situation described above. In this situation, the temperature T 3 , 2 is derived from the sum 

of the gravitino decay rates involving the coupling constants A' 132 , A' 211 , A 223 , A 311 , Am, 

A211, A233 an d A323. 

We see in Fig.f|(f ) (dashed line) that the temperature T 3 , 2 fulfills the condition |H] but not 



14 In fact, gravitino masses larger than ~ 30TeV are needed so that the condition [12] can 
be fulfilled. This can be observed on Fig.|5| (dashed line) in which have been represented 
the values of the gravitino and superpartner masses for which kT 3 , 2 > OAMeV. Those val- 
ues have been derived within the optimistic scenario described above. In this optimistic 
scenario where several fl p coupling constants are simultaneously present, the gravitino 
masses corresponding to kT 3 , 2 > OAMeV are smaller than in the cases of a single domi- 
nant fl p coupling constant, since as we have already seen T 3 , 2 is typically proportional 
to the fl p coupling constants and decreases with the gravitino mass. It is also interes- 
ting to note that in this optimistic scenario, for a value of the ratio g{T d ) l l A / ' giT^,^) 1 ^ 2 
(see Eq.(|9])) of order 50, the condition [L^ is fulfilled for LSP gravitino masses of 0{TeV) 
so that the constraint m ~ 0(TeV) due to the hierarchy problem can be respected. 
However, for g{kT^/ 2 w 5 10~ 3 MeV) ~ 5 (see Eq. flTTl) ), corresponding to the value of 
the SU(3) C x SU(2) L x U(1) Y model, g(T d )^ 4 / 'g(kT 3/2 w 5 lO^MeV) 1 / 2 w 50 leads to 
g(Td) ~ 10 8 which is not realistic within the present models. 

In summary, if the gravitino mass, and thus the scalar superpartner masses, do not 
exceed the TeV scale by at least one order of magnitude, the gravitino decay temperature 
T 3 / 2 fulfills well the constraint of Eq.(^) but not of Eq . ([12]) , within the scenario of an 
unstable LSP gravitino decaying via A, A' or A" interactions. Hence, this scenario does 
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not appear to be a natural solution to the cosmological gravitino problem. 

5 Conclusion 

Along with the existence of the cosmic microwave background, big-bang nucleosynthesis 
is one of the most important predictions of the big-bang cosmology. Furthermore, if one 
assumes that the light nuclei (atomic number less than 7) have effectively been produced 
through the big-bang nucleosynthesis, one finds that the theoretical predictions on the 



abundances of these light nuclei are in good agreements with the observational data p7 
Now, as we have seen above, if one believes that the light elements have been synthesized 
through the standard big-bang nucleosynthesis, the gravitinos must have decayed before 
the nucleosynthesis epoch. Therefore, since we have found that this cannot happen in 
the scenario characterized by an unstable LSP gravitino having a decay channel which 
involves trilinear ftp interactions, this scenario does not seem to provide a realistic solution 
to the large relic abundance of the gravitino. It must be specified that this conclusion has 
been obtained by assuming that the gravitino and thus the scalar superpartner masses 
do not exceed the TeV scale so that supersymmetry can solve the so-called hierarchy 
problem, which is one of the strongest motivations for the existence of supersymmetry. 
We also mention that this final conclusion is based on the present constraints on the 
trilinear fl p couplings JTTJ] and on the assumption that the ratio g(Td) 1//4 / 'g(Ts/2) 1 ^ 2 takes 
reasonable values. 



The authors of |TJ| have found that the lifetime of the partial two-body gravitino decay 
G — > 7 + v, initiated by the bilinear ftp interactions with a coupling constant yielding 
contributions to the neutrino mass matrix which match the observed atmospheric neutrino 
anomaly, reads as: 



, 26 ,'lGeV\» 



r 3 / 2 (G — 7 + v) « 8.3 10 2b sec, y\x>) 

\ m 3 / 2 J 

where the neutrino-photino mixing matrix element was set at the value appropriate to a 
photino with a mass of 80GeV. If one assumes tentatively that the relation in Eq.([H|) con- 
tinues to hold in order of magnitude for the heavier gauginos case, then the condition that 
LSP gravitinos decay radiatively before the nucleosynthesis epoch, namely r 3 / 2 ~ 10 2 sec, 

requires the lower mass bound m 3 / 2 ~ 2 10 8 GeV which conflicts with the constraint from 
the hierarchy problem. On the other hand, the case of light LSP gravitinos, say with 
mass m 3 / 2 < lGeV, leads within the bilinear ftp option to very long lived gravitinos 
which could constitute a valid candidate for the dark matter 1113 . 



Within the ftp trilinear interactions option, by invoking the bounds on coupling constants 
rather than the values which match the observations involving the neutrinos mass ma- 
trix, we obtain values of the gravitino lifetime which are significantly shorter than those 



obtained within the bilinear interactions option by Takayama and Yamaguchi |LJ] . As a 
function of the relevant ratio of squared masses of the gravitino and exchanged sfermions, 
z = my 2 /rh 2 , the predicted partial gravitino lifetime, associated to the decay into a 
given flavor configuration of three light fermions final state, decreases monotonically in 
the interval < z < 1 as, 

/ 1 \ /]TpV\ 3 
r V 2(G -> fififk) « 10 9 - 10 11 (j-)( sec, (19) 

K ^ijk J V m 3/2 / 
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where A stands for any one of the trilinear coupling constants. For light gravitinos of mass 
m 3/2 = OiGeV) and Ay* = 0(1), we see from Eq. ([T9|) that the gravitino lifetime exceeds 
by a few order of magnitudes the age of the universe to — 3.2 10 sec. 

We have found that the present bounds on the trilinear fl p couplings allow easily the 
gravitinos to decay before the present epoch. In consequence, the scenario of an unstable 
LSP gravitino having a decay channel which involves trilinear fl p interactions could pos- 
sibly solve the cosmological gravitino problem in the context of an inflationary model 



28| . The reason is that within such a framework, the gravitino would not have necessa- 
rily to decay before the nucleosynthesis epoch, since in inflationary models small gravitino 
number densities can be predicted (for small enough gravitino number densities, the gra- 
vitino decays do not upset the nucleosynthesis). Indeed, if the universe has experienced 
inflationary expansion, the primordial gravitino abundance would have been completely 
diluted by the exponential expansion. Then, after the universe would have been reheated, 
gravitinos would have been regenerated by scattering processes off the thermal radiation, 
resulting in a secondary gravitino number density proportional to the reheating tempe- 
rature, which is a parameter of the model. 
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A Definitions of the functions involved in the gravitino decay 
rates 

The mass functions Q(m), E(?77.i, 777,2), A(?77.i, 777,2) and E (777,1, 777.2) used in Section ||] are 
defined in the following, 

fi( m ) = — f 60m 6 - 162m 4 m 2 /2 + U0m 2 my 2 - 37m 6 /2 

+ 12(m 2 - m§ /2 ) 3 (5-^- - 1) log(l - ^)) , (A.l) 
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5(777.1,777.2) = — ( 84m 1 m 2 m 3 / 2 + 7m 3 , 2 + tt7 1 (— 12m 1 + 24m 2 + 42m 1 m 3 / 2 — 40m 



3/2; 



+m 2 (-i2m 2 y + 24m^ + 42m 2 m 2 /2 - 40m£ /2 ) + 12m 2 ( m ^ m ^ — H!l + 4 m * 

m 3/2 
2 

— 4m 2 m 3//2 — hm\m\i 2 + 2m 3//2 + 2m 2 m 2 ) log(l ^-) + 12m 2 

( — | + 4m 2 — Am 2 m 2 i 2 — hrn\m\/ 2 + 2m 3 / 2 + 2m 2 m 2 ) 

m 3/2 

logd-^f)), (A.2) 

2 2 

A(mi, m 2 ) = -m\m\{ x 2 2 + m 2 + m 2 - m| /2 ), (A.3) 

^ m 3/2 

2 o 2 

S(mi, m 2 ) = A(mi, m 2 ) log(l - ^) log(l + ^ - ^), (A.4) 

777,1 m 2 m 2 

m3/2 being the gravitino mass. 

B Formulas for the spin summed amplitudes of the gravitino 
decay reactions 

The full analytical result of the scared amplitude summed over the spins for the gravitino 

decay reaction G -^ ViCjCk is the sum of the following scared amplitude and interference 
terms 1 , 

1 A 2 

I ji/f |2 x ijk / 2 2 1 2 \f 2 2 2 \ 

1 a| = 3 M 2 (m 2 -m 2 p ( m 3/2 ~ ™,fc + ™,JKfc ~ ™ ej ~ ™ e J 

2 \2 



( m 3/2 + m ?fc 



mr 



"3/2 -T- " *jfc "»„<; 2» m n 

-2 m jfc , (B.1J 



4m 3/2 

1 A 2 

I jit 1 2 L ijk I 1 2 i 2 \ / 2 2 2- 

' 6 ' = 3 MHm 2 - ml ) 2 ^ 3/2 ~ ifc "V^ ~ m * ~ m ^ 



lr The projective sums for a massive spin-3/2 field are require d in ord er to calculate the spin averaged gravitino decay 
amplitudes. These projective sums are available in the literature |12| , [28[ BST . 
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where M a , M^ and M 3 / 2 correspond to the amplitudes of the three contributions repre- 
sented in Fig.Q and marked by the same letters a, b and c. In Eq. (|B.l|) to Eq. (|B.6|) , we 
have used the same notations as in Section ^| and the quantities m^-, rrij k and m^ are 
defined by, 



m- 



(jpM+piej))'' 
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m jk = (p(ej)+P(e*)) 

ml = (P^)+P(e k )) 2 , (B.7) 

p{ u i)i Pi e j) an d p(efc) being respectively the quadri-momenta of the particles z/j, e^ and e k . 



^ A 



The formula for the spin summed amplitude of the gravitino decay process G -^ Vi€je k , 
which is given in Eq. (p.l| )-Eq.( p.6| ), is equal to that of the charge conjugated process 

G — > ViCjCk. 

The full analytical result of the scared amplitude summed over the spins for the gra- 

A' _ ~ A' \ 

vitino decay reaction G -^ Vidjd k {G -^ CiUjd k ) is equal to that of the reaction G -^ 
ViZjZk given above, simply by changing \j k to A^- fe , adding an extra multiplicative color 
factor of N c = 3 and replacing the masses as follows: rno iL — > m„ iL (m e - iL ), mg. L — > 
m d ]L ( m u 3L ), m ShR -> mg kR (m SkR ), m Vi -»■ m Ui (m e J, m e . -> m dj . (m Ui ) and m ek -> 
7^4 ( m d fc )- The definitions of the quantities my, m^fc and m^ must also be changed into, 



«4 = ipiy%) + p{dj)f 

m% = (p(dj) + p{d k )f 

ml = (p(^)+p(4)) 2 , (B.8) 

for the spin summed amplitude of the reaction G -^ Uidjd k , and into, 

m 2 ij = (P(ei) + P(uj)) 2 
m % = (p(uj) + p(d k )) 2 
ml = Wi) + P{d k ))\ (B.9) 

for the spin summed amplitude of the reaction G -^ eiUjd k . The same remark as above 

holds for the charge conjugated process. 

The full analytical result of the scared amplitude summed over the spins for the gra- 

~ K k 
vitino decay reaction G -^ Uidjdk is the sum of the following scared amplitude and 

interference terms, 
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(B.12) 



where M a and M& are associated to the amplitudes of the two contributions represented in 
Fig.||] and marked by the same letters a, b and c. In Eq.( [B.ld|) to Eq.( |B.12|) , the notations 
of Section |3] have been used and the quantities m^-, rrijk and m^ are defined as follows, 



(B.13) 



K 


= {p{ui)+p{dj)Y 


™% 


= (p{dj)+P(dk)) 2 


m'ik 


= (p(Ui)+ P (d k )) 2 . 



The above remark on the charge conjugated process still holds. 



2i) 
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